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Finite size analysis of the pseudo specific heat in SU(2) gauge theory 

J. Engels and T. Scheideler'* 

'^Physics Department, University of Bielefeld, P.O. Box 100131, 33501 Bielefeld, Germany 



We investigate the pseudo specific heat of SU(2) gauge theory near the crossover point on 4* to 16^ lattices. 
Several different methods are used to determine the specific heat. The curious finite size dependence of the 
peak maximum is explained from the interplay of the crossover phenomenon with the deconfinement transition 
occurring due to the finite extension of the lattice. In this context we calculate the modulus of the lattice average 
of the Polyakov loop on symmetric lattices and compare it to the prediction from a random walk model. 
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1. INTRODUCTION 

The pseudo specific heat Cy of SU{2) gauge 
theory was already investigated in the begin- 
ning of Monte Carlo lattice studies. It is known 
to have a peak near (3 — w 2.2, in the 

crossover region between strong and weak cou- 
pling behaviour. A first finite size analysis by 
Brower et al. ^ on 4* — 10^ lattices revealed a 
strange dependence of the peak on the volume 
V = {N^aY of the lattice. Here, a is the lat- 
tice spacing and N„ the number of points in each 
direction. The location of the peak shifts with in- 
creasing volume from smaller /3— values to larger 
and then to smaller ones again; the peak max- 
imum decreases with increasing volume. Such 
a behaviour is unknown for any ordinary phase 
transition. The nature and origin of the peak 
remained therefore unclear, though a connection 
to the nearby endpoint of the first order critical 
line in the (/3, /3— adjoint)-plane was proposed by 
Bhanot and Creutz Recently, new calcula- 
tions for this extended SU{2) model were per- 
formed [^[^]. Likewise, a new study of Cy with 
higher statistics and also on larger lattices, utiliz- 
ing the analysis techniques now available, seems 
appropriate. In addition, since symmetric lattices 
are used to simulate zero temperature physics, it 
is important to estimate remaining finite temper- 
ature effects which may be seen in the pseudo 
specific heat. 
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2. METHODS TO CALCULATE Cy 

We use the standard Wilson action for SU{2) 



where 



(2) 



is the plaquette or energy and [/^,y(a;) is the 
plaquette link operator. The sum extends over 
all independent forward plaquettes. There are 
Np = 6N^ such plaquettes. We denote the lattice 
average of the plaquettes by P 



The speudo specific heat is then defined by 



Cy = 



d{P) ^ _g2<m 
d{l/0) ' d(3 
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There are three methods to determine Cy : 

i) one measures the plaquette expectation values 
(P) as a function of (3 and calculates the numer- 
ical derivative at (3m = P + A(3/2 from 

CviPM) = -^i{P)i[3 + A(3) - (F}(/3)) ; (5) 

ii) one measures the variance of the plaquettes, 
which is proportional to Cy 



Cy = P'NpUP^) - (Pf) ; 



(6) 
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iii) one calculates the sum of plaquette-plaquette 
correlations 



{P?) 



(7) 



As it should be, all three methods are in complete 
consistence with each other. The most straight- 
forward way is, of course, to calculate the variance 
of (P). The density of states method (DSM) may 
then be used to interpolate between the points. 
The variance becomes, however, definitely too 
small, if the plaquettes are measured during and 
not after each update, because that leads to local 
correlations among the plaquettes. 
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Figure 1. The pseudo specific heat Cy calcu- 
lated from the variance ( squares), the numer- 
ical derivative (circles) and the DSM interpola- 
tion, calculated on a 6'' lattice. 



In Fig. 1^ we show for comparison the results 
from methods i) and ii) for a G'* lattice. Here, on 
the average 90-120 thousand measurements were 
made every fifth update at each value. An 
update consisted of one heatbath and two overre- 
laxation steps. 

We have also investigated the plaquette cor- 
relations. We find in general a rapid fall with 
R = x' — X, the correlation length is of order 
1. The plaquettes Pf_iu{x) and Pfi'u'(x') may be 
in parallel or orthogonal planes. At the peak 



(/3 K, 2.23) we find that the total contribution of 
the orthogonal correlations is about 30% higher 
than that of the parallel correlations, whereas far 
away from the peak, at /3 = 2.70, the contribu- 
tions are essentially equal. 

3. FINITE SIZE DEPENDENCE OF Cv 

In Fig. H we compare the results for Cy from 
lattices with N„ = 4, 6, 8, 12 and 16. The general 
behaviour already found in [|l[ is fully confirmed. 
On the other hand, we see that there is no fur- 
ther finite size dependence in the peak region, 
if Ncr > 8. This suggests, that the finite size 
dependence of the smaller lattices is related to a 
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Figure 2. The pseudo specific heat Cy vs. (3 on 
lattices. Part (a) shows the results for — 4: 
and 6 from the DSM interpolation, for No- = 8, 12 
and 16 the measured points were connected by 
straight lines. Part (b) shows the peak region for 
N„ = 8, 12 and 16 in detail. 
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different phenomenon. Indeed, the critical point 
for the Nt = 4 finite temperature deconfine- 
ment transition is at fic = 2.30, very close to the 
crossover peak positon. Since we are using pe- 
riodic boundary conditions for all directions, the 
approach to the critical point corresponding to 
Na will influence the plaquette expectation val- 
ues. To check this, we have calculated the lattice 
average L of the Polyakov loop L{x) 

1 1 ^' 

^ = E ^(^) ; ^(^) ^ 2 n f^-.* • (8) 

" X t = l 

As can be seen from Fig. H, the expectation 
value of the modulus of L is not zero on sym- 
metric lattices, not even in the strong coupling 
limit, i.e. we have finite temperature effects also 
on symmetric lattices. Well below j3c{N„) the 
quantity (|L|) is a constant. With increasing j3 
it starts to increase already before the transition 
point. It is obvious, that due to the nearby tran- 
sition points the crossover peaks of the 4^ and 6^ 
lattices are stronger distorted than those of the 
larger lattices, where only the right shoulders of 
the peaks are slightly influenced. We may find the 
A^cr— dependence of sA (3 — Q from a simple 
random walk model, where the modulus of a 
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Figure 3. The expectation value of the modu- 
lus of L for N„ = 4,6,8, 12 and 16 vs. (3. The 
broken vertical lines show the locations of the cor- 
responding finite temperature phase transitions. 



sum over A^^^ equal random variables is propor- 

3/2 

tional to No- ■ Indeed, we obtain from a simula- 
tion at /3 = the following relation for SU{2) 
(|i|)/3=o - 0.400 • A^-3/2 

The values, which we calculated at /3 = 2.0 

are shown in Fig. |^ together with the results for 
/3 = 0. At P = 2.0 all lattices apart from the 4^* 
lattice have already reached the strong coupling 
value. Deviations from this value indicate then 
the onset of finite temperature effects. 

Our final conclusion is, that the crossover peak 
is not the result of an ordinary phase transition. 
For large lattices the peak is at /3co = 2.23(2), its 
height is Cv,co = 1.685(10). 
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Figure 4. The logarithm of (|L|) vs. InA^o- on 
4-^ - le'^ lattices at /3 = (sohd fine), at /3 = 2.0 
( circles ) and from an estimate from the {\L\) 
value at Ncr = 1 (dashed line). 
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